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ABSTRACT 



An experimental investigation was performed to determine, 
in situ, the change in mean particle size across the exhaust 
nozzle of a small metallized solid-propellant rocket motor. 
Light scattering profiles were recorded at both the exhaust 
and the entrance of the nozzle. The experimental method 
used provided excellent results within the exhaust. However, 
combustion light at the wavelength of the transmitted light 
hampered light scattering measurements within the motor. 
Particle size measurements were consistent with the sizes 
found in the collected exhuast products . 
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I. 



INTRODUCTION 



The importance of aluminum oxide (AZ^O^) particle size 
produced in solid-propellant rocket motors has been recog- 
nized since aluminum was first introduced as a fuel additive. 
Aluminum added to a propellant will increase the thrust 
produced and aid in suppressing high-frequency combustion 
instability. Although the performance improves with aluminum 
addition, two phase flow losses adversely affect specific 
impulse efficiency. Specific impulse efficiency is the product 
of several component efficiencies, i.e., chemical-kinetic, 
combustion, two-phase flow, two-dimensional flow, submergence 
and throat erosion. In aluminized propellants, the two-phase 
flow loss is one of the more significant losses. Two-phase 
flow losses result from particle-gas thermal and velocity lags. 
The latter strongly depends on particle size. Parametric 
studies involving small motors, with throat diameters less 
than one inch, at AFRPL [Ref. 1] and CSD [Ref. 2] have 
revealed that a 1 micron change in particle diameter results 
in a change of 1% in Isp. The effective use of particle 
damping for suppression of combustion pressure oscillations 
is also strongly dependent upon the particle size distribution. 

One of several analytical computer programs used to pre- 
dict performance of a solid-propellant rocket motor is the 
SPP (Solid Performance Program) [Refs. 3,4]. As the above 
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discussion indicates, it is of critical importance to accurately 
predict particle size. This is especially true for particle 
damping with minimum smoke propellants which typically con- 
tain less than 2% metal by weight. The 1975 version of the 
SPP [Ref. 3] contains a particle sizing empirical model 
developed by Cohen. This model was based on correlating 
collected exhaust samples. After a critical review, this model 
was considered inadequate. A new model was therefore required 
for a new SPP version, titled Improved SPP [Ref. 4]. Inorder 
to meet SPP accuracy goals, this model had to predict particle 
size within ±10%. Unfortunately, there are no adequate 
theoretical models relating particle size to propellant and 
motor parameters. This necessitates the use of empirical 
methods for correlating particle size. 

There are two methods that can be employed to develop an 
empirical relation. One is to correlate particle size data to 
propellant and motor variables, as was done by Cohen. The 
other method is to calculate particle size based on a critical 
Weber number for a maximum stable drop size [Ref. 5]. After 
considering the relative merits of each method, Hermsen [Ref. 

5] devised an empirical model based on collected particle size 
data for the improved SPP. This model was correlated to a much 
broader data base than was previously available to Cohen. 

The particle diameter correlated was , the mass-weighted 
diameter, given by: 
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43 



y n.d 4 / y n. d 3 

1 l' L 11 



Smith [Ref. 6] has shown that was the most appropriate 

average to use in predicting two-phase flow losses. 

Hermsen's empirical relation [Ref. 5] is presented 
below : 



D 



43 



3.6404 



. 29 32 . . - . 000 8163£Pt 
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where : 

= mass-weighted average diameter, 

= nozzle throat diameter, in. 

E, = A £ 2 °3 concentration in chamber, g-mol/10 0 g 

P = chamber pressure, psia 

t = average residence time, msec 



The predicted is now applied to the expression below in 

order to calculate the two phase flow loss [Ref. 7]. 



ETATP 



C 3 ( 5C 4 d Vp' 15 £ 




) 



where : 



ETATP 



D 



P 



impulse loss effect, % 

mole fraction of condensed species, moles 
of condensed species per 100 grams of 
mixture 

in microns 
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p 



average motor chamber pressure, psia 



e 



nozzle expansion ratio at ignition 
conditions 



D 



t 



nozzle throat diameter at ignitiion, 
inches 



e 



Naperian base 2.71828 



C( ) 



coefficients that are a function of D 
and D 



t 



P 



The two-phase flow loss is applied with the other aforementioned 
efficiencies to a theoretical Isp. The Isp loss from boundary 
layer effects is then subtracted from the corrected theoreti- 
cal Isp to obtain delivered Isp. 

Collecting exhaust products is feasible only for small 
motors and the techniques used often result in variation in 
the measured size. Validation of the SPP model must be based 
on data collected within the motor environment. Light 
scattering measurements are well-suited for this application. 
These measurements are more easily made if the particles are 
large compared with the wave length of the transmitted light 
and the transmittance is greater than 90%. The former allows 
the accurate application of Fraunhofer diffraction, whereas the 
latter specification is necessary in order to satisfy single 
scattering requirements . 

The Naval Postgraduate School has conducted a series of 
investigations to determine the applicability of light scatter- 
ing measurements to the solid rocket motor exhaust products 
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[Refs. 8,9,10,11]. This investigation was directed at modi- 
fying the apparatus and techniques developed in the previous 
efforts in order to more accurately determine particle size. 
Validation of the apparatus was also a goal. The SPP particle 
prediction model can then be verified by a comparison of the 
predicted size to the experimental data. 
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II. THEORETICAL BACKGROUND 



The most general theory of light scattering was developed 
by Mie and is presented by Van de Hulst [Ref. 12]. Mie's 
theory, although not restrictive in its assumptions in regard 
to particle size and refractive index, is however very com- 
plex to utilize for data reduction. Dobbins, et al. [Ref. 13] 
and Powell, et al . [Ref. 14] derived similar expressions for 
intensity ratios at two forward scattering angles based upon 
an Upper Limit Distribution Function (ULDF) . and Fraunhofer 
diffraction. The ULDF was initially proposed by Mugele and 
Evans in Reference 15. Although the expressions developed by 
Dobbins and Powell are extremely useful in relating intensity 
ratios to D ^ > the volume-surface mean diameter, these models 
are still somewhat cumbersome to use for data reduction 
purposes . 

Particle sizes found in solid-propellant rocket motors 
are often bi-modal . Particles in the larger mode are the most 
important for two-phase flow loss calculations and are gener- 
ally much greater in diameter than the wavelength of inci- 
dent light. Scattering by these large particles is described 
by Fraunhofer diffraction. Hodkinson [Ref. 16] and others 
[Refs. 13,14] have demonstrated that diffraction theory 
(within the forward lobe) can be used to approximate the Mie 
theory. Buchele in Reference 17 summarized the experimental 
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techniques for determining the particle size by measurements 
of dif fractively scattered light. Buchele presents a function 
which closely approximates the curves given in References 13 
and 14. 



1(0) = EXP [-(.57 a 0 ) 2 ] 



where : 



1(0) = ratio of intensity of scattered light at 

some angle (0) to the intensity of 
scattered light at theta equal to zero 
degrees 

a = ttD/X is the particle size parameter for 

diameter D _2 and wavelength of light 
lambda (A) 7 



The above expression was used in the current study for data 
reduction to obtain Unfortunately, is the mean 

diameter used in the SPP. Therefore, roust be related to 

either by analyzing exhuast samples or by approximate 
analytical means. 

A final important point is that the above expression is 
valid for a <_ 3 (the center lobe) . For the expected ranges 
of theta and D^' this proves to be a very mild restriction. 
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III. EXPERIMENTAL APPARATUS 



A. ROCKET MOTOR 

The rocket motor used in this investigation was the same 
as the long motor described in Reference 11. The motor com- 
ponents are presented in a photograph in Figure 3.1 and the 
installed motor is shown in Figure 3.2. A schematic is also 
presented in Figure 3.3. The propellant grain dimensions re- 
mained the same as in the previous effort [Ref. 11]. The 
grain was two inches in diameter, one inch in length and had 
a web thickness of .725 inches. The grain was cylindrically 
perforated and allowed to burn also on the aft end. The 
igniter was composed of BKNO^ • Ignitiion of the igniter was 
achieved by heating a resistance wire with a 12 VDC source. 

Some of the exhaust products were collected in an eight 
inch diameter stainless steel exhaust tube. The products 
were collected from both the side walls of the tube and the 
end cap. Another particle sample was collected from the walls 
of the combustor, near the nozzle entrance. Both samples were 

V 

then cleaned and mounted for observation using a scanning 
electron microscope. 

B. LIGHT SCATTERING APPARATUS 

The apparatus shown in Figures 3.4 and 3.5 was basically 
the same as used in Reference 11 with some exceptions. The 
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light source used for the nozzle exhaust was a 8-miliwatt, 
helium-neon laser. The nozzle entrance also used a 8-milliwatt 
helium-neon laser. The lasers were mounted on two parallel 
optics benches, such that one beam could pass through the ex- 
haust and the other beam could pass through the motor cavity 
adjacent to the nozzle entrance. A beam expander/collimator 
was used with the helium-neon laser in the exhaust. Each beam 
was intercepted by a stop placed directly in front of the 
narrow pass filter. The narrow pass filter was used to eliminate 
extraneous light. The stop and narrow pass filter were placed 
30.5 centimeters from the motor centerline for the exhaust 
beam and 11.5 centimeters from the motor centerline for the 
motor cavity beam. The exhaust beam narrow pass filter had 
a diameter of 5.08 centimeters and the motor cavity narrow 
pass filter had a diameter of 2.3 centimeters. Both beams 
passed through the centers of these filters. In the present 
arrangement of the apparatus this resulted in a maximum angle 
for collecting scattered light of 4.8 degrees in the exhaust 
and 6.0 degrees in the motor cavity. The placement of the 
optics closer to the motor centerline allowed for a greater 
maximum angle to be recorded than in the original design. 

Directly following the narrow pass filters were condensing 
lenses. Both paths used a lens with a diameter of 5.08 
centimeters. These lenses had a focal length of 50 centimeters 
and were utilized to focus scattered light onto a linear array 
of photodiodes. The array consisted of 1024 silicon photodiodes 
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with .025 millimeter spacing. The apparatus in Reference 11 
positioned the lenses such that only angles from 0 to .04 
radians could be observed. In actuality, however (as a conse- 
quence of diffraction around the beam stop and beam spreading) , 
the range of useful angles where only scattered light occurred 
was from .02 radians to .04 radians. In the current arrange- 
ment the diode array was placed 1.0 centimeters below the main 
beam and the collecting tube diameter was enlarged. This 
allowed a range of angles from .02 to .07 radians to be 
observed . 

C. DATA ACQUISITION AND DATA REDUCTION 

The data acquisition system is presented in detail in 
Reference 10. The system controller was an HP 9836S computer 
and the HP 6942A multiprogrammer provided A/D conversion and 
storage. The data acquisition system was programmed to take 
data at a specified point in the motor firing by inputing a 
trigger pressure and a time delay after the trigger pressure 
was reached. The number of diode scans for the exhaust and 
motor cavity were eight and four respectively (limited by 
the current data acquisition system memory size) . 

Two methods were used in References 10 and 11 for data 
reduction, one of which was altered for this study. The modi- 
fied method pertained to the iterative graphics technique. 

In References 10 and 11 the experimental data were normalized 
and compared to a theoretical profile that was also normalized 
to some assumed scattered light intensity at zero degrees 
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A "best-fit" of the experimental profile 



for a specified D^. 
to the theoretical one yielded The current investigation 

determined* that this method resulted in ambiguity in deter- 
mining D ^2 as a consequence of assuming a centerline intensity 
for the scattering profile. The method used in the present 
study was one in which the experimental data were not altered. 

A theoretical curve was generated for any specified value of 
D ^2 based upon a selected initial angle (0^) and the correspond- 
ing intensity at that angle (1^) from the unaltered experimental 
data. The approximate intensity profile for a polydispersion 
given by Reference 17 is : 

I 2 = I 1 EXP-{ (0 2 - 6^) ( .57 Dr/A) 2 } 

Specification of a series of values for 0 2 results in 
values for I 2 and, therefore, an I vs 0 theoretical profile 
that passes through the experimental data at I, , 0]_. If 
this curve does not fit the experimental data, then other values 
of D ^2 can be chosen and a new theoretical curve generated. 

A "best-fit" method is subsequently used to determine D^* 

In summary, the method used for obtaining the mean particle 
size is based upon fitting an approximate (although quite 
accurate within specific angle restrictions) theoretical 
profile for a polydispersion to the unaltered profile obtained 
experimentally over a small range of scattering angles. The 
modified data reduction program is presented in the Appendix. 
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The second method was similar to the first and also used 



the above expression solved for D^ 2 : 

Dj 2 = -Ln(I 2 /I 1 ) (A/.57TT) 2 /(02 - 0^) 

The procedure used was to first select a starting angle, 0^, 
along with a 0 2 /0^ ratio. A value of I 2 at © 2 was obtained 
from the experimental data. This was repeated with increasing 
values of © 2 and then again with different values of 0 2 /0^. 

The results were plotted graphically as particle size vs 
scattering angle for each angle ratio. A nearly constant 
value of D 22 vs 0 indicated a good correlation. 
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TABLE I 




L ^ J 


r 

Description 


Copper 

Nozzle 


Copper 

Nozzle 


Outside Diameter 
(d, inch) 


2.125 


2.125 


Length (In, inch) 


1.25 


1.25 


Throat Diameter (dt, inch) 


.28 


. 23 


Slope (converging) 
Angle (y, degree) 


45 


45 
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Figure 3.1. Motor Components 
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Figure 3.2. Installation of Rocket Motor 
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Figure 3.3. Schematic of the Motor 
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Figure 3.4. Light Scattering Apparatus 
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Figure 3.5. Schematic Diagram of Light Scattering Apparatus 



TABLE II 



LIGHT SOURCE SPECIFICATION 
A. Helium-Neon Laser (Exhaust) 



1. Manufacturer 


Spec tra-Phy sics 


2 . Model : 


147 


3 . Type : 


He-Ne Class IIIB 


4. Output Power: 


8 mWatt 


5. Beam Diameter: 


. 92 mm 


6. Beam Divergence: 


.87 nrad. 


B. Helium-Neon Laser (Motor 


Cavity) 


1. Manufacturer: 


Uni-Phase 


2. Model: 


1305P 


3 . Type : 


He-Ne Class IIIB 


4. Output Power: 


5 mWatt 


5. Beam Diameter: 


. 81 mm 


6. Beam Divergence: 


1.00 mrad. 
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IV. RESULT AND DISCUSSION 



A. INTRODUCTION 

The purpose of this study was to improve and validate 
the accuracy of the apparatus used for particle sizing in 
the rocket motor environment. Also, the experimental data 
will provide a basis to verify the particle size prediction 
model of the SPP. 

The apparatus had been calibrated by Harris [Ref. 10] 
and Kertadidjaja [Ref. 11] using glass and polystyrene spheres 
in a suspension of water. Kertadidjaja [Ref. 11] also 
attempted to demonstrate the applicability of the apparatus 
to the combustion environment. This investigation continued 
the earlier studies with several refinements. Three propel- 
lant types were used in this investigation. All three utilized 
a GAP binder with ammonium perchorate (AP) oxidizer. Each had 
varied aluminum content, with the percentages present being 
0%, 2%, and 4.8%. The size of the aluminum powder was 20 
microns. The major ingredient compositions for the 4.8% 
aluminized propellant were 14,. 7% GAP, 45.7% 200 pm AP, 24.6% 

25 pm AP, 8.5% TEGIN and 4.8% 20 pm Al . The 0% and 2% 
aluminized propellants kept the same solids loading. 

B. SYSTEM CALIBRATION 

As previously mentioned, the apparatus was calibrated by 
both Harris [Ref. 10] and Kertadidjaja [Ref. 11]. The 
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apparatus was re-calibrated in this study as a result of the 
modifications that were made. The standard calibration proce- 
dure described in detail by Harris [Ref. 10] remained the 
same. The laser beam passed through a glass container holding 
a suspension of polystyrene spheres and water and was initially 
positioned on the number one diode. The linear diode array 
was then moved vertically down, with the aid of a dial indica- 
tor, such that the number one diode was one centimeter below 
the main beam. Since the diode array was 2.56 centimeters 
long and was 50 centimeters from the focusing lens, this 
arrangement allowed for a range of measuring angles from .02 
to .07 radians. Therefore, the scattered light striking the 
diodes would not be affected by the main transmitted beam or 
the light diffracted by the beam stop as in the previous two 
studies. This arrangement was used both for the exhaust beam 
and the motor cavity beam. 

Two polydispersion samples of polystyrene spheres, with 
D 32 ' s ^.8 microns and 10.2 microns, were used for the cali- 

bration. These samples were chosen because the D^ 2 was known 
and the expected °f the actual aluminum oxide particles 

during a motor firing was between 1 and 15 microns. The 
results are presented in Figures 4.1, 4.2, 4.3 and 4.4. These 
results were encouraging in regards to validating the modifi- 
cations and the method for data reduction. Calibrations were 
also conducted within the motor cavity with the laser beam 
and scattered light passing through the motor windows. Figures 
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4.1 and 4.2 are representative of the results obtained. This 
indicated that the motor windows did not adversely affect 
scattering profile and, therefore, bias the results. 

An initial point of concern was the narrow range of measur 
ing angles which were employed in the present apparatus. 

Since large particles diffract most of the incident light at 
small angles and very small particles diffract most of the 
incident light at large angles, the range and minimum angle 
employed could not have resulted in sufficient sensitivity to 
particle size for the expected range of particle sizes. Also 
shown in Figure 4.1 are two additional theoretical profiles 
for D ^2 °f 10 and 8.3 microns. It was apparent from the sensi 
tivity of the slope change with that the mean particle 

size could be determined with an uncertainty less than ±.5 
microns . 

The "upsweeps" in the curves at larger angles using the 
two angle method were the result of the increase in slope of 
the scattering profile due to the effect of truncation by 
the apparatus. These parts of the curves remain only to 
provide a reference to the angle ratios. 

Another informative point -in apparatus design is given 
by Hirleman [Ref. 18] in regards to placement of the focusing 
lens. Hirleman notes that truncation of the larger scatter- 
ing anbles, which is a function of the placement of focusing 
lens, will cause a biasing against the small particles. The 
expression given in this discussion is: 
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(ttcI/7 . 66 A) D 



Z 

c 



where : 

Z = critical distance in front of the focusing 
lens where significant energy diffracted 
by particles less than diameter D is 
truncated by the focusing lens 

A = wavelength of the incident light 

For a given Z , all particles that are. greater than D in size 
will have 84 percent of their diffracted energy pass through 
the focusing lens. In the exhaust configuration of the present 
apparatus was equal to 30.5 centimeters, d was 5.08 centi- 
meters and A was .6328 microns. This yields a D of 9 microns. 
Thus, according to this criterion, some of the scattered light 
from particles below 9 microns is truncated and, therefore, 
results in a bias towards larger particles. This type of 
biasing was small in the present apparatus down to particle 
sizes of approximately 3 microns (the smallest particles in 
the 4.8 micron polydispersion) as evidenced by the good 

calibration results is shown in Figure 4.3. However, the 
apparatus could be improved by either moving the lens 
closer to the exhaust (at the risk of apparatus damage from 
the exhaust jet) or by using a larger diameter lens. 

Transmittance effects were also investigated in conjunc- 
tion with the calibration. As previously mentioned, trans- 
mittances above 90% are required to meet single scattering 
requirements. Calibrations were conducted using polystyrene 
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spheres with a of 10.2 microns. Scattering data were 

recorded for transmittances of 85%, 70%, 60%, 50%, and 30%. 

The results are presented in Figures 4.5 through 4.14, and 
summarized in Table III. Both the curve fit and the two- 
angle methods resulted in very good agreement with the actual 
mean particle size with transmittances as low as 50%. Even a 
transmittance of 30% did not severely affect the measurement. 
These results indicate that particle sizing methods based upon 
single scattering theory can effectively be used in a multiple 
scattering environment with transmittances significantly less 
than 90%. 

A test was then conducted to determine the transmittance 
in the actual motor exhaust using the 2% aluminized propellant. 
The procedure was the same as the one used in Reference 11. 

The transmitted beam was directed onto the diode array (with 
the appropriate filters to reduce the intensity) and the beam 
profile was recorded. The motor was then fired and the profile 
was again recorded. The results are presented in Figure 4.15. 
The results indicated a transmittance greater than 90%, with 
some beam spreading. There was also no apparent beam shift. 
Therefore, transmittance loss was determined not to be a 
debilitating factor in the motor exhaust jet. 

The calibration curves have an apparent limit of approxi- 
mately .04 radians for which scattered light can be recorded 
(for example, see Figure 4.3). However, it was stated that 
the apparatus was designed to record angles up to .07 radians. 
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This is due to the water-glass interface having an effective 
index of refraction of 1.39. Thus, the scattering angles are 
corrected for this high index of refraction, resulting in a 
corrected maximum angle of approximately .04 radians. 

C. NOZZLE EXHAUST MEASUREMENTS 

One test was conducted with a nonmetallized propellant. 

The results, presented in Figure 4.16, show no scattering. 
Therefore, the scattering observed by Kertadidjaja [Ref. 11] 
at smaller angles for a similar test was probably the result 
of beam broadening and diffraction around the beam stop. 

Three tests were conducted with both the 4.8% and 2% 
aluminized propellants. The motor firing data and a summary 
of the particle size measurements are presented, respectively, 
in Tables IV and V. From light scattering measurements each 
propellant showed a similar result, one test with a of 

approximately 5 microns and one with a larger value of 
(7-12 microns). This test-to-test variation probably resulted 
from the very short burn times and the corresponding lack of 
a well-defined steady-state condition. These results are 
presented in Figures 4.17 through 4.28. Typical S.E.M. photo- 
graphs are shown in Figure 4.29. The light scattering results 
were quite consistent with those obtained from the S.E.M. 
evaluation. Although the sample sizes from which the S.E.M. 
evaluations were based were relatively small (see Table VI) , 
it still provided a good indication as to the validity of 
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the light scattering measurements. It should be noted that 
^> 22 ' an< ^ to an even greater extent D^, is strongly dominated 
by the larger particles. One or two large particles can 
significantly change 2 • Future S.E.M. evaluations should 
utilize much larger sample sizes. 

The SPP equation for (which is primarily a function of 

throat diameter) predicted much smaller values of than 

observed in this limited set of data. The SPP model, however, 
is based on particulate data from motors with nozzle diameters 
greater than 1 inch. 

D. NOZZLE ENTRANCE MEASUREMENTS 

Light scattering measurements through the motor at the 
nozzle entrance were unsuccessful as a result of combustion 
light existing at the same frequency as the laser light. This 
was the case for both the helium-neon laser (A = .6328 ym) 
and the argon laser (A = .488 ym) . 

The analyses of the S.E.M. data for the nozzle entrance 
indicated a much larger and as compared to the nozzle 

exit measurements. Therefore, a significant reduction in size 
appears to have occurred to the larger particles across the 

it 

exhaust nozzle. 
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TABLE III 



EFFECT OF TRANSMITTANCE 
ON MEASURED PARTICLE SIZE 
(D ^2 from Manufacturers Data Equals 10.2) 



Transmittance 

(%) 



Measured D^ 
(microns) 



85 



9.8 



70 



9.6 



60 



9.5 



50 



9.4 



30 



9.3 
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TABLE IV 



SUMMARY OF MOTOR FIRINGS 



Date of 
Test 


Wt of 
Alum 
(%) 


Press 

Pc 

(psig) 


Press 
Max P 
(psig) 


Burn 

Time 

(sec) 


Residence 

Time 

(msec) 


D f t 

(in 


Oct 31, 85 


0 


350 


366 


3.8 


- 


.28 


Oct 22 , 85 


4 . 8 


312 


328 


2.9 


10.9 


. 28 


Oct 23, 85 


4.8 


328 


359 


2.3 


11.3 


.28 


Oct 18, 85 


4.8 


333 


340 


2.8 


11.3 


. 28 


Nov 19, 85 


4.8 


450 


466 


. 95 


14.0 


. 25 


Oct 30, 85 


2.0 


296 


296 


2. 8 


10.4 


. 28 


Nov 26, 85 


2 . 0 


408 


467 


2.0 


14 . 4 


. 28 


Nov 25, 85 


2.0 


566 


575 


1.8 


16.7 


. 25 



Note: Pressure colume labeled P c indicates pressure 

at which data was taken. 
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TABLE V 



SUMMARY OF EXPERIMENTAL RESULTS 



Wt of 
Alum 
(%) 


Press 

Pc 

(psig) 


d 32 

Light Set. 
(psig) 


d 32 

S.E.M. 

Ext./Mtr. 


d 43 

S.E.M. 

Ext./Mtr. 


D 43 

SPP Model 


4.8 


328 


5.6 


5.6/18.2 


7.3/24.7 


2.3 


4.8 


333 


5.5 


- / - 


- / - 


2.3 


4.8 


450 


12.0 


13.6/18.9 


19.7/23.9 


2.5 


2.0 


296 


7.0 


12.2/14.4 


18.6/19.6 


2.3 


2.0 


408 


9.3 


- / - 


- / - 


2.3 


2.0 


566 


4.6 


4.9/10.8 


6.0/13.2 


2.5 



Mote: (-) indicates sample was not taken for this run. 

Diameters are in microns. 



Ti 
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TABLE VI 



S.E.M. PARTICLE DISTRIBUTION 



Oct 23 Ext. 
Oct. 23 Mtr. 
Oct. 30 Ext. 
Oct. 30 Mtr. 
Nov. 19 Ext. 
Nov. 19 Mtr. 
Nov. 25 Ext. 
Nov. 25 Mtr. 



12 3 4 

53 22 16 11 

50 21 22 13 

45 42 11 13 

24 15 3 6 

64 35 31 7 

73 29 9 5 

44 39 22 3 

37 30 11 1 



Size in mien 

5 6 7 8 

7 4 2 1 

9 6 8 4 

21 2 2 

28 3 7 13 

6 11 1 

12 4 1 1 

14 4 1 1 

15 8 3 8 



9 10 12 15 

1 

12 11 
1 1 
6 3 

2 1 

12 11 
1 

3 4 3 



19 20 25 30 

1 2 
1 1 

12 1 
1 1 

1 2 1 

1 
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CURVE FIT RESULTS 
INTENSITY vs . THETR 
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Figure 4.1 Calibration With 10.2 Micron Polystyrene Spheres 
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Figure 4.2 Calibration With 10.2 Micron Polystyrene Spheres 
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Figure 4.3 Calibration With 4.8 Micron Polystyrene Spheres 
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Figure 4.4 Calibration With 4.8 Micron Polystyrene Spheres 
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Figure 4.5 Curve Fit Method, 85% Transmitted Light 
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Figure 4.6 Two Angle Method, 85% Transmitted Light 
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Figure 4.7 Curve Fit Method, 70% Transmitted Light 
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Figure 4.8 Two Angle Method, 70% Transmitted Light 
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Figure 4.9 Curve Fit Method, 60% Transmitted Light 



THO-RNGLE METHOD 




OJ — * 

(suojoloi) DZIS 



49 



THETR , (rad) 

Figure 4,10 Two Angle Method, 60% Transmitted Light 
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Figure 4.11 Curve Fit Method, 50% Transmitted Light 
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Figure 4.12 Two Angle Method, 50% Transmitted Light 
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Figure 4.13 Curve Fit Method, 30% Transmitted Light 



t?9£0 * 




OJ — • 

( suojo LUJ) 3ZIS 



53 



THETR (rad) 

Figure 4.14 Two Angle Method, 30% Transmitted Light 
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DIODE NUMBER 

Figure 4.15 2% Aluminum Propellant, Less Than 90% Attenuation 
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Figure 4.16 Non-Metalized Propellant, 22 Oct., Exhaust 
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Figure 4.17 4.8% Aluminum Propellant, 18 Oct., Exhaust 
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Figure 4.18 4.8% Aluminum Propellant, 18 Oct., Exhaust 
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Figure 4.19 4.8% Aluminum Propellant, 23 Oct., Exhaust 
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Figure 4.20 4.8% Aluminum Propellant, 23 Oct., Exhaust 
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Figure 4.21 2% Aluminum Propellant, 30 Oct., Exhaust 
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THETR (rad) 

Figure 4.22 2% Aluminum Propellant, 30 Oct., Exhaust 
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THETR (rad) 

Figure 4.23 4.8% Aluminum Propellant, 19 Nov., Exhaust 



THO-RNGLE METHOD 




CM 



T3 

rd 

E_ 




(SUOJOLUU) 



3ZIS 






63 



Figure 4.24 4.8% Aluminum Propellant, 19 Nov., Exhaust 
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THETR (rad) 

Figure 4.25 2% Aluminum Propellant, 25 Nov., Exhaust 
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Figure 4.26 2% Aluminum Propellant, 25 Nov., Exhaust 
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Figure 4.27 2% Aluminum Propellant, 26 Nov., Exhaust 
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Figure 4.28 2% Aluminum Propellant, 26 Nov., Exhaust 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The results of this study have shown that the measure- 
ment of D ^2 in rocket motor exhausts can be made using quite 
simple light scattering methods. Apparatus calibration 
results and S.E.M. evaluations of the collected motor exhaust 
products support this conclusion. 

The SPP equation model underestimated D^ for this inves- 
tigation. This indicated that the model predictions cannot 
be extrapolated to small motors with throat diameters less 
than one inch. The light scattering method, however, can 
be used to validate the SPP model for the larger motors upon 
which it was based. 

In order to improve the signal to noise ratio in regards 
to the scattered light intensity, it is recommended that the 
laser beam be moved closer to the exit plane of the nozzle, 
where the particle density is greater. The laser beam for 
this study was 7.3 centimeters aft of the nozzle exit. To 
resolve the combustion light problem it may be necessary to use 
a different laser, such as a CO^ that has a different wavelength 
than present in the combustion process. Also, other binder 
systems should be tested to determine if they significantly 
alter the light spectrum. 
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APPENDIX 



PROGRAM LISTINGS 



10 

20 

30 

40 

50 

60 

70 

80 



i »x mxxxxx xxmxm xxxx*»***x**xxx***m*»*»x**xx»x xxxmxxmxxxxx 



iimmiiimimiii 
i **»*»»xx*x»x*«»xx*x 
!**»x**x»xxxxxxxxx*« 
i xxxxx) ix*«xx*»xm»x 
! xxxxxmxxmxxxx** 
lllSOMIHIHmiHli 



ROC 

PLOTS RAU DATA 
AVERAGES 
FILTERS 

DETERMINES MEAN DIAMETER 
BY INTERACTIVE GRAPHICS 
AMD THE TWO-ANGLE METHOD 



i»***»*xmxx*m*r* 

**x*xx*xmxx**x**»» 

»X »»3HXXX»XXX»m»X X 

**x*x*x****x»xxx*xx* 

*XJX#»XXmX»X*>X»X* 

XX**IIX»XXXX»XXXX»*X» 

*xxxxixx»xxxxx»*x«x* 



90 i*xxxxx»*x»xxxxx*x»* Robert Kelly Harris 84 xxxx»»*xxxx*xxx*x*xx 

100 '»<»x*»»»»*x»*mxxxx Revised, John S Rosa 85 ****»»x*»*»*x*xmxx 

110 !»***mmx*x**»*»»****M*»»****m»*»»»**mm»»»»»*»*x*m»x»»»* 

120 OPTION BASE 1 

130 COM /Hrdaaus/ Av2(1024) 

131 COM A2U024) PUFFER 

140 COM /Gauss/ TK 1024), L 

150 CCM /Max/ M7,M5, Xt,Yt,Xr*,Yn,Xx 

160 COH /Readata/ B,P,H,g3$l20],Q4»l20),Zz4I201,Yl(81?2) BUFFER ,Y2(81?2) BHFFE 
R 

170 COM /Two/ A*K1024),H,H1,F 

18? DIM Scans (B ) , X 1 1 034 ) 

190 COM /Plots/ P14I201 ,P2$I20] ,P3tl201,P4tl201 

200 INTEGER Graf ( 1 : 12480) BUFFER 

210 Choose: PRINT CHRt(12) 



280 

230 

ON’ 

240 

258 

260 

270 

280 

2 ? 0 * 

300 

310 

320 

330 

340 

350 

360 

370 

380 



Old=0 

PRINT * TO LOOK AT ANY DATA FILE THE PROGRAM NEEDS SOME STARTING INFORHATI 
PRINT ’’ 

PRINT " TO ACCOUNT FOR THE CHANGE IN WAVELENGTH IN THE MEDIUM’ 

PRINT ” 

PRINT ’ ENTER THE INDEX OF REFRACTION OF THE MEDIUM’ 

PRINT ’’ 

PRINT ’ WATER = 1,33" 

PRINT’ AIR = 1,0" 

PRINT’ ESTIMATE OF EXHAUST = 1.1 ’ 

BEEP 

INPUT ’ THIS VALUE ADJUSTS THE COMPARISON CURVE TO THE MEDIUM’, Ml 

PRINT CHR$(12) 

PRINT * TO ACCOUNT FOR REFRACTION OF LIGHT AT BOUNDARIES BF.TUEF.N’ 

PRINT 

PRINT ’ THE MEDIUM AND AIR YOU MUST ENTER THE INDEX OF REFRACTION ■ 
PRINT ’’ 
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PRINT ■ OF THE COMBINATION OF THE MEDIUM AND IIS BOUNDARY' 

PR T NT " 

PRINT ' THIS VALUE IS APPLIED DIRECTLY TO THE DATA' 

PRINT " 

PRINT ' AIR=t .0 PLEXIGLASS 6 WATER = 1.39' 

PRINT " 

PRINT 'ESTIMATE OF MOTOR CAVITY COMBUSTION PRODUCTS 6 WINDOW = 1,22* 

PRINT " 

PRINT ' ESTIMATE OF ROCKET EXHAUST * 1.1 OR 1.0' 

INPUT M 

PRINT CHRK 12) 

INPUT ' ENTER LASER WAVELENGTH (F?Np=. 6328 ,Ar=. 483) \L 

INPUT 'ENTER THE DESIRED PLOTTING INTERVAL OF DIODE ARRAY (2,4,6)\H 
F=500 !nn FOCUL LENGTH OF OBJECTIVE LFNS 

X<1 

Tt(1)=,92 

Diod=.025 im DIODE SPACING MAY BF .03 

FOR 1=1 TO 1023 

X ( 1+1 )=I+1 1CREATF AN ARRAY OF DTODF NUMBERS 

Tl(I+l)=ATN((I*Diod)/(F*M)MTl(l) 'CREATE AN ARRAY OF THETA 

NEXT I 

PRINT CUR! M2) 

CALL Displavl(Old) ! JUST PRINTS INITIAL REMARKS ON CRT 

IF Old 1HEN GAIL RenpwlAvl 00) 'THESE TWO LINES ALLOW FUR 

IF Old THEN GOTO Gauss ! REVIEWING REDUCED DATA FILES 

FF.EP 

INPUT ' INPUT TWO FILENAMES NOW (NO-PART ,PART)",Q3$,Q41 

INPUT * ENTER '\' FOR 8 SCANS FOR 4 SCANS* ,P 

7;1-':IN1FRNAL,4,1‘ ! USES LEFT DISK DRIVE FOR RAW DATA FILES 



390 
400 

410 
420 
430 
440 
450 
460 
470 
480 
491) 

500 
510 

520 

521 

522 
530 

549 

550 
560 
570 
580 
590 
600 
6T0 
620 
630 
640 
650 

669 CALL Rpadata 

670 CALL Shift 

680 BEEP 

6 7 0 Serpen: ' 

700 PRINT CMRi (12) 

710 P 1 l-'VOLTAGE' 

720 P2<=”D10DF NUMBER" 

730 P 31 -"VOLTAGE us. DIODE' 

740 P4$-"RAW DATA’ 

759 M5=1024 

751 Mfl^O. 

760 Xt -64 

770 Y t= . 04 

780 Xm=4 

790 Y«=5 

791 X*=0 

800 CALL Plot(Y2(»)) 

910 GSTORE GrafU) 

820 f/LL Dataplot(B,Yl(*),H) 
830 PRINT ' 

840 PRINT "' 

850 PRINT ' 

860 INPUT J 



'CLEARS SCREEN 
! STRINGS FOR PLOTS 



! SET UP VALUES FOR PLOTS 



! DRAWS AXES ,FTC... 

! STORES GRAPHICS D ISLAY JUST MADE 
1PLOT NO-PARTICLES DATA 
NO-PARTICLES PLOTTED' 

HOW MANY SCANS SEEM TO BE GOOD 7" 
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B O PRINT ■* 

I PRINT 1 WHICH SCANS ARE GOOD .ie. . .1 ,2, 4,5,7, INCLUDE LAST CO 



MMA' 


ILAST COMMA IS REQUIRED OR YOU HAK TO HIT CONT.. TWICE 


890 


INPUT Scans!*) 




900 


PRINT CHR4U2) 




910 


PRINT USING */////////' 




920 


PRINT * 


AVERAGING THE SELECTED SCANS' 


930 


GRAPHICS OFF 




940 


CALL Average! J, Scans!*), Yl(*),Avl(*)) 


950 


PRINT USING *9* 




960 


GLOAD Graf!*) ILOADS GRAPHICS ARRAY RATHER THAN WASTE TIME RE-DRAWING 


970 


GRAPHICS ON 




980 


CALL Result!Awl(*),X(*),H) 


! PLOTS 1024 ELEMENT ARRAYS 


990 


PRINT USING '//' 




1000 


PRINT * 


Average Intensity No-Particles' 


1010 


ON KEY 0 LAPEL 'AVERAGE!' 


GOTO Screen 


1021 


ON KEY 1 LABEL 'PLOT-PARTICLES* GOTO 1060 


1030 


PRINT ' 


PRESS KEY 0 0 TO RE-AVERACE OR * 1 TO CONTINUE' 


1040 


Standby: ! 


MANUAL CALLS THIS INTERRUPT DRIVEN PROGRAMMING 


1050 


GOTO Standby 


! LOOPS, WAITING FOR USER TO DECIDE 


1060 


PRINT USING 




1070 


OFF KEY 1 


! HELPS AVOID CONFUSION BY CLEARING THAT POX 


1080 


CALL Datapl«t(B,Y2(*),H) 


! PLOT PARTICLES DATA 


1090 


PRINT USING '/' 




1100 


PRINT ' 


PARTICLE DATA PLOTTED* 


1110 


PRINT ' FOR A HARD COPY OF THIS RAW DATA PRESS KEY * 6" 


1120 


PRINT * ' 




1130 


PRINT ' 


OR TO CONTINUE PRESS KEY 1 1" 


1140 


ON KEY 6 LABEL 'HARD i RAW 


' GOTO Rau 


1150 


ON KEY 1 LABEL 'CONTINUE " 


GOTO Select 


1160 


GOTO Standby 




1170 


Raw: CALL Plot (Y2( ») , 1 ) 


! THE ONE (1) IS AN OPTIONAL PARAMETER 


1180 


CALL Datapl#t(B,Y2(*),H) ! WHICH IS USED TO GET HARD COPIES 


1190 


Select:! 




1200 


OFF KEY 6 


“ 


1210 


OFF KFY 1 




1220 


PRINT CHRi! 12 ) 




1230 


PRINT USING '////' 




1240 


PRINT ' 


HOU MANY SCANS SFEM TO BE GOOD ?' 


1250 


INPUT J 




1260 


PRINT ' 


WHICH SCANS ARE GOOD .ie... 1,2, 4, 5, 7, INCLUDE 


LAST COMA* 




1270 


INPUT Scans!*) 




1280 


PRINT CHR$(12) 




1290 


PRINT USING '/////////' 




1300 


GRAPHICS OFF 




1310 


PRINT ' 


AVERAGING THE SELECTED SCANS' 


1320 


CALL A*erage(J,Scans(*),Y2(*),Av2(»)) 


1330 


PRINT USING "8* 




1340 


P4$=*AVERAGED SCANS' 
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1350 CALL Plot(A»2(*)) 

1360 GRAPHICS ON 

1370 GSTORE Graft*) 

1380 CALL Result<Av2(«),X(*),H) 

1390 PRINT USING *//' 

1400 PRINT * Average Scattered Intensity* 

1410 OH KEY 2 LABEL ’AVERAGES* GOTO 1080 



1420 ON KEY 3 LABEL ’SUBTRACT* GOTO 1450 

1430 PRINT * HIT KEY * 2 TO RE-AVERAGE OR KEY 0 3 TO 

CONTINUE* 

1440 GOTO Standby 
1450 PRINT USING *8* 

1460 OFF KEY 3 

1470 MAT Avl= Av2-Avl (SUBTRACTS NO-PARTICLES FROM PARTICLES 

1430 GLOAD Graf!*) 

1490 CALL Result(Avl(*),X(*),H) 

1500 PRINT USING "///* 

15t0 PRINT * Plot of the Difference Between Particles and N 



O-Partulps" 



1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 



PRINT * KEY * 6 FOR HARD COPY* 

PRINT * KEY * 1 10 CONTINUE " 

ON KEY 6 LABEL * HARD AVERAGE* GOTO 1570 
ON KEY 1 LABEL * FILTER * GOTO 1600 
GOTO Standby 



CALL PlottAvl (*), 1 ) 

CALL Result(Avl(*>,X(*),H) 
PLOTTER IS 3 f “INTERNAL" 
PRINT USING *//* 

PRINT * 

°RINT " 

PRINT * 

PRINT * 

INPUT * 

IF F i 1=0 THEN Gauss 



ENTER THE NUHBFR OF TIMES YOU WISH TO APPLY" 
THE DIGITAL FILTER FCR SMOOTHING" 
EXAMPLF******** 10" 

TAKES ABOUT 1.5 MINUTES* 

ZERO IS ALSO ACCEPTABLE", Fil 



P4t=*FlLIERED DATA* 

CALL FiIteMAv1(*),Fil) 

PRINT CHRtt 12) 

CALL PlottAvK*) ) 

IF Old THEN MAT X= T1 
CALL ResulttAolt*) ,X(*),H) 

PRINT * Plot of the Difference Between Particles and N 



O-Particlps" 



1740 PRINT 
TEP* 

1750 
1760 
1770 
1780 
1790 
1800 
1810 



PRINT * 

PRINT * 

PRINT * 

PRINT * 

PRINT * 

ON KEY 8 LAPEL 
ON KEY 9 LABEL 



AFTER APPLICATION OF A DIGITAL FIL 

for a HARD COPY * 
PREPARE the PLOTTER and* 
PRFSS KEY I 9* 

OR* 

PRESS KEY * 9 TO NORMALIZE* 



'HARD FILTERED* GOTO 1830 
1 NORMALIZE * GOTO Gauss 
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! The 1 is for hard copy 



1820 

1830 

1840 

1850 

I860 

1870 

1880 

1870 

1891 

1892 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2920 
2030 
'2031 
2040 
2050 
2060 
2070 
2080 
2090 
2190 
2120 
21 30 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 



GOTO Standby 
C8LI Plot(Avl(»),l) 

CALL Result! AvU*), XU), H) 
OFF KEY 8 
Gauss: OFF KEY 2 
OFF KEY 8 
OK KEY 1 LABEL ’ 



! KEY 1 OPTION IS STILL VALID 



! SET UP VALUES 
'FOR PLOTS 



FILTER ■ GOTO 1600 ! ALLOWS ONE TO SMOOTH OLD DATA 
! WHEN REVIFUING REDUCED FILES 

01d=l ! THIS ALLOWS FILTERING EVEN AFTFR ENTERING GRAPHICS NODE. THE 

i DIODE COORDINATES AREN'T NEEDED ANYMORE »SFE LINES NEAR CALL FILTER 
PRINT CWW12) 

1 Av2 Array is Normalized BUT AvlU) is Still Saved For Re-work if Needed 
MAT Av2= Avl ! TO BFGIN WITH THE ARRAY IS ASSUMED TO 9E NORMALIZED 
! AND THE USER CHANGES THIS WITH INTERACTIVE GRAPHICS 
i IN THE SUBROUTINE "Cotopare" 

PM=’FILTFRED VOLTAGE 

P29="THETA (rad)" (STRINGS FOR PLOTS 

P?t="INIENSITY vs, THETA" 

P 4 t- "CURVE FIT RESULTS* 

M5=,070 
Xt= .005 
Yt=.l 
Xm=2 
Ym=2 
Xx=, 02 

CALL P)ot(AvlU)) 

CALL Result(AvlU) ,T1(*),H) 

GSTORE Graf!*) 

CALL Comp are! Avl(*),H,H, HI, Graf!*)) 

GRAPHICS OFF 

ON KEY 4 l ABEL "OTHER ARRAY" GOTO Choose 

ON KEY 5 1 ABEL " TWO-ANCLE" GOTO Buchele 

ON KEY 6 LABEL "STORE DATA" GOTO Storedata 

ON KEY 7 LABEL * QUIT * GOTO Quit 

CALL Hen eel 
GOTO Standby 
Buchele: 1 

CALL Two angle 
CALL Menu el 
GOTO Standby 
Storedata: 1 
CALL Store 
CALL Menuel 
GOTO Standby 



(PLOT OF NORMALIZED INTENSITY vs. 
! SAVFS THE SCREEN IMAGE 



THETA 



! PRINTS OPTIONS ON THE SCREEN 



Quit: 

END 



i 



SUB Averane(J, Scans!*) ,YU),Av(*)) (AVERAGES SELECTED SCANS 
HAT Av= (9) 'INITIALIZES THE ARRAY LOCAL TO THIS ROUTINE 

FOR 1=1 TO J STEP 1 

K=(Scans(I)-l)*1024+l ! THIS COUNTER IS THE BEGINNING 

L=Scans(I)*1024 ! AND THIS ONE THE END OF BLOCKS 

FOR 11 =K TO L (OF 1924 IN THF. OVERALL DATA 
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2320 

2330 


12=11 -XU ■ THIS COUNTER IS ALWAYS PF TWEEN 

Av(I2)=A<i(I2W(It) ! 1 AND 1024 


2340 


NEXT 11 




2350 


NEXT I 




2360 


HAT Av= A«/(J) 




2370 


SUBENO 




2380 


SUP Datao]ot(B,Y(*),H) 


!B IS 4 OR 8 (THE NUMBER OF SCANS) 


2390 


COM /Max/ M7,H5,Xt,Yt,X«,Y(i ) Xx 




2400 


l DIR 0 




24t0 


LORG 4 




2420 


FOR J=1 TO P 


!EACH SCAN 


2430 


MOVE 500+J*50,H7- .05 


! SEE NEXI LINE 


2440 


LAPEL Jj ! HELPS KEEP TRACK OF FACH SCAN AS IT APPFARS 


2450 


LINE TYPE 1 




2460 


K=(J-1)*10?.4U 


! BEGINNING OF EACH SCAN AND 


2470 


L-IX1024 


•THE END WITHIN THE TOTAL BLOCK 


2480 


HOVE 1 ,Y(K) 


! MOVE TO THE FIRST POINT 


2490 


FOR T=K TO L STEP H 




2500 


X-I-K+1 


! THIS GIVES 1 TO 1024 FOR ABSCISSA 


2510 


PLOT X, Y( I) 




2520 


NEXT I 




2530 


NEXT J 




2540 


PENUP 




2511 


M8=0 




2550 


SUBEND 




2560 


SUP Resu]t(Iav(*),D(*),H) 


! THE AVERAGE INTENSITY IS PLOTTED 


2570 


CLIP ON 




2580 


MOVE DU) ,Iav(1 ) 




2590 


FOR 1=1 TO 1024 STEP H 




2600 


DRAW DU),lav(l) 




2610 


NEXT I 




2620 


PENUP 




2630 


STIPEND 




2631 


Soli: ! 




2640 


SUB Plot(YU), OPTIONAL Hard) 


• IF THE OPTIONAl.(Hard)IS PFCEIVFD 


2650 


COM /Max/ H7,M5.Xt ,Yt ,X*,Y«,Xx 


'THE FIGURE GOES IO THE PLOTTER 


2660 


COM /Plots/ Pit 1201 ,P2tl201 ,P3il20I ,P4t(201 


2670 


M7=MAX(Y(*>)#1.1 


! A SCALING VARIABI E 


2680 


M6=.I*M7 


! ANOTHER SCALING VARIABLE 


2700 


C1NIT 


'GINIT IS JUST GOOD PRACTICE SI) YOU KNOW 


2710 




! WHERE YOU ARE BEGINNING 


2720 


SELECT NPAR 


! THIS DFIFCTS IF THE HARD COPY TS DESIRED 


2730 


CASE 1 


1 1T COULD BE DONE WITH IF THEN LOGIC 


2740 


PLOTTER IS 3 /INTERNAL" 


! BUT IS PRESEN IED FOR FAMILIARIZATION 


2750 


CASE 2 


ISINCE IT IS HCRE POWERFUL IN COMPLEX 


2760 


PLOTTER IS 705/HPGL* 


(SITUATIONS 


2770 


END SELFCT 




2790 GRAPHICS ON 
2790 CSIZE 5,. 6 
2800 DEG 
281 0 LDIR 0 


! TO BE ABLE TO SEE THE PLOT 
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2820 LO»G 5 
2R30 HOVE 75 ,9*5 

2840 LABEL USING 'K';P4* 

2850 HOVE 75,99 
2860 LABEL P3$ 

2870 HOVE 75,20 
2880 LABEL P2$ 

2890 LOIR 90 
2930 HOVE 18,60 
2910 LABEL PI* 

2920 VIEWPORT 30,125,30,95 
2930 FRAHE 

2940 WINDOW Xx,H5,-H6,H7 

2950 AXES Xt,M6,Xx ,0,Xn,I ,5 

2960 LOIR 0 

2970 CSIZE 3,. 5 

2980 LORG 8 

2990 CLIP OFF 

3000 FOR I=-H6 TO H7 STEP H7/10 

3019 H8VE X»,I 

3020 LAPEL USING ’# ,HD. DDDD’ } I 

3030 NEXT I 

3040 CSIZE 3,. 6 
3050 LDIR 90 

3069 LORG 8 

3070 FOP I=Xx TO H5 STEP Xt 

3080 HOVE T ,-H6 

3090 LAPEL USING ’t,K*;I 

3100 NEXT I 

3110 PENUP 

3120 SUBEND 

3130 SUB Conparo(Av1(*),H,H, HI, INTEGER Graft*)) 
3140 CQH /C«ss/ T1(«),L 

3150 COH /Hrdqaus/ Av2(») 

3160 COH /Max/ H7,H5,Xt,Yt,Xn,Yn,Xx 

3201 Pool: ! 

3202 PRINT * INPUT REFERENCE ANGLE * 

3293 INPUT Tt 

3205 FOR N=t TO 1024 STFP 1 

3206 IF T1(N))=Tt THEN 3208 

3207 NEXT N 

3208 lt=0 

3299 FOR I=N-5 TD N+5 

3210 It=It+Awl (I) 

3211 NEXT I 

3212 Al=It/l 1 

3220 OFF XEY (GETS 

3239 ON KEY 3 LABEL "HENUE* GOTO 4160 

3240 PLOTTER IS 3, 'INTERNAL* ! IN CASE 



! SCREEN UNITS FOR MARGINS 
! DRA’JS A BOX 



INUHBER THF Y AXIS 



INUHBER THE X AXIS 



RID OF ALL LABELS ON KEYS 
A HARD COPY WAS JUST HADE 
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3250 Hard=0 ! (0) SO ONF. DOESN'T EXIT TOO SOON 

32S0 Centerline=l ! THE INITIAL NORMALIZING VALUE 

32.70 D=10. ! INITIAL PATICLE MFAN DIAHETER IN MICRONS 

3280 Change: CSIZE 4 ,.6 

3290 PRINT USING *//////* 

3300 PRINT * STAND BY FOR CURVE" 

3310 HOVE 0,1 

3320 CLIP ON 

3331 DIM T(1:1024) 

3340 FOP 1=1 TO 1024 STEP H 

3350 Tl I)=AI*EXP(-D A 2*( (Tl< I) A 2-Tt A 2)*(( (Ml* , 57«PI )/L) A 2> ) ) 

3370 DRAW T1(I),T(I) 

3375 IF (Tl (I)*( (PI*D*M1 )/D) >3,0 THEN 3440 

3380 NEXT I 

3440 PENUP 

3450 IF Hard THEN 3550 IEXITS ROUTINE IF A HARD COPY HAS JUST MADE 



3469 PRINT USING "/////////" 

3470 PRINT " 

SIZE" 

3480 PRINT * 

3490 PRINT " " 

3500 PRINT ’ 

3510 PRINT " 

3520 PRINT * * 

3530 PRINT " 

3540 PRINT USING *//////" 

3541 OFF KEY 0 

3550 ON KEY 6 l ABEL "HARD COPY" GOTO 

3560 ON KEY 9 LARFL "NORMALIZE" GOTO 

3570 ON KNOB .05 GOTO Poise 

3580 Wait: GOTO Wait 
3590 Pulse: PRINT USING "f»* 

3690 Strnnl4="Mean Dianeter " 

3610 Strnq2$=* Microns* 

3620 C»unt=KNOPX 

3630 D=DR0UND(D+Count/15,2) 

3640 GLOAD Graft*) 

3650 LDIR 0 

366 0 LOIG 8 

3670 MOVE M5,.8*M7 

3680 LABEL Strrgl$6VAL*(D)iStrng2$ 

3690 GOTO Change 

3701 Hardqauss: ! 

3710 PRINT " 

3720 PRINT " 

3730 PRINT " 

3740 PAUSE 

3750 Hard=l 

3760 CALL Plot (Av2<») ,1 ) 



USE THE KNOB TO VARY THE PARTICLE 

OR HIT KEY 4 6 FOR HARD COPY* 

KEY \ 9 ALLOWS YOU 10- 
SELECT NEW REF ANGLE* 

HIT KEY i 3 TO GET OUT" 



Hardqauss 

Pool 

! ( .05) IS TIME INTERVAL IN WHICH 
! PULSES FROM THE KNOB APE 
! COUNTFD AMD THIS NUMPEP IS 
! USED BY THE INTERRACTIVE 
! GRAPHICS TO VARY THE PARTICLE 
! SIZE AND PLOT THE ASSOCIATED 
! GAUSSIAN APPROXIMATIONS OF 
! SCATTERING PROFILES 



PREPARE the PLOTTER* 
PRESS CONTINUE for" 
HARD COPY* 

! SO THAT SUPEXIT OCCURS AFTER HARD COPY 
! 1**HA«D COPY 
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377(1 


CALI Result(Av2<»),Tl(«),H) 


3780 


PRINT USING ’ 8 ' 


■ 


3790 


LORG 8 




3300 


HOVE H5,.8»M7 




3810 


LDIR 0 




3820 


LABEL StrnqmVALKD)4Strng2i 


3800 


GOTO Change 




4160 Subexit: ! 




4170 


SUBEND 




4180 


SUB Hpnuel 




4190 


PRINT USING 


"( 5 * 


4200 


PRINl " 


YOU CAN RE-NORMAL I7.E' 


4210 


PRINT x" 


4220 


PRINT n 


YOU CAN RE-AVERAGE (New Data Onlv)" 


4230 


PRINT ■xxxxxxxxxxxxxxxxxxxmxxxxxxxxxxxxmxxxxxxxxxxxxxxxxxHxxxx* 


4240 


PRINT ’ 


PRESS KEY ♦ 4 TO LOOK AT OTHER DATA* 


4250 


PRINT • 


OLD OR NEW 1 


4260 


PRINT ’ 


MOTOR / EXHAUST’ 


427| 


PRINT ■xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxixxxxxxxxxxxxxxx* 


4280 


PRINT ■ 


PRESS KEY 1 5 * 


4270 


PRINT * 


FOP THE ' TUO-ANGI E ' METHOD" 


4300 


PRINT ' 


OF PARTICl E SUING" 


4310 


PRINT "xxxmxxxxxxxxxxxxmxxm#mxxxx#mx*xyxx*xxmxx#xxxxxx* ,, 


4320 


PRINT " 


TO STORE THE REDUCED DATA PRFSS KEY \ 6 


4330 


PRINT "xxxxxxxxxxxxxxxxxxxxmxxxxxxxxxxxxxxxxxxxxxrxxxxxxxxxxxxxxx-* 


4340 


PRINT * 


TO END THIS SESSION PRESS KEY 1 7" 


4350 


SUBEND 




4360 


SUB DisplavKOld) 


4378 


PRINT * 


TO REDUCE NEW DATA PRESS KEY \ 1* 


4380 


PRINT ** 




43?0 


PRINT ' 


TO REVIEW PREVIOUSLY REDUCED DATA PRESS KEY » 2* 


4400 


PRINT ’ 1 




4410 


ON KEY 1 IABPL 1 


NEW DATA’ GOTO New 


4428 


OK KEY 2 LAPEL ’ 


OLD DATA" GOTO Review 


4430 Wait: GOTO Wait 




4448 


New: PRINT CHR$(12 


) 


4450 


PRINT * 


PUT THE DISK IN THE LEFT DRIVE IF IT IS NOT ALREADY’ 


4460 


PRINT ” 




4470 


PRINT * 


ENTER THE NAMES OF THE FIRST AND SECOND FILES.’ 


4480 


PRINT "• 




4490 


PRINT 1 


EACH FILE HAS DATA FROM BOTH DIODF ARRAYS.’ 


4500 


PRINT ■■ 




4510 


PRINT ’ 


FIRST IS NO-PARTICLES 6 NEXT IS PARTICLES Dlt,D2$’ 


4520 


PRINT ” 




4530 


Old=0 




4548 


SUPEXIT 




4550 Review: 01d=l 


(THIS VARIABLE IS PASSED TO THE MAIN PROGRAM TO INDICATE 


4560 




•THAT THE DATA TO PE READ IS ONE (1) SCAN AND THAT 


4570 




(THE AVERAGING ROUTINES WILL NOT BE USED 
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4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 

4720 



SUP-END 
SUB Readata 

CPM /Readata/ B,P ,H,Q3*120 1 ,04M201 ,Zz*[20 ] ,Y1 (*) BUFFER, Y200 B'OFFFR 



Xyz=l 
P=8 

R 1=65536 

IF P=2 THEN Xy z=40?7 
IF P=2 THEN R 1=32768 
IF P=2 THEN B=4 
PRINT " * 

PRINT ” 

ASSIGN RDi=kl TO Q3WZzl 
ASSIGN PDisk2 TO Q4$6 Zj$ 
ASSIGN RB'iff 1 TO BUFFER Y1 (») 
ASSIGN 8Buff2 TO BUFFER Y200 



! FILE POINTFR VAR [API E 
! NUHBFR OF SCANS IN MOTOR DA (A 
(NUMBER OF BYTES OF MOTOR OAF A 
! RECORD 4 WHERF 4 SCANS DATA BEGINS 
! NUMBER OF BYTES OF 4 SCAN DATA 
(NUMBER OF SCANS 



READING DATA FROM FILE ON DISK" 



(OPEN I/O PATHS 



4730 CONTROL BDi=kl ,5;Xyz 

4740 TRANSFER PDiskl TO PBuffl [COUNT R1 

4750 WAIT FOR EOT PDiskl 

4760 CONTROL PDisk2,5;Xyz 

4770 

4780 TRANSFER POisk2 TO ?Buff2;C0UNT R1 

4790 WAIT FOR EOT PDisk2 

4800 ASSIGN PDiskl TO » 

4810 ASSIGN PDisk? TO # 

4820 ASSIGN ?P«ffl TO * 

4830 ASSIGN EBuff2 TO * 

4840 SUBEND 
4850 SUB Plot2(F,D(*),X) 

4860 DIN Titlel*l25],ThetaM201 

4870 GINIT 

4880 IF X=1 THEN PLOTTER IS 705,'HPGL’ 

4890 DEG 

4900 GRAPHICS ON 

4910 VIEWPORT 35,125,35,85 

4920 Max=10*INT((MAX<D<»))+10)/10) 

4930 WINDOW 0(1,2), DIE, 2), 0, Mai 

4940 G=(D<E,2)-D(1 ,2»/(E-l)#4 

4950 F=( INK E/4) -1 ) »4 

4960 IF F=0 THEN F=4 

4970 CLIP D(4,2)-G,D(F,2)+2«G,0,Max 

4980 GRID G,5,D(4,2)-G,0 

4990 CLIP OFF 

5000 LORG 8 

5010 LDIR 90 

5020 CSIZE 4, ,5 

5030 FOR 1=4 TO E STFP 4 (PUTS NUMBERS ON X AXIS 



(DEGREES For LABEL DIRECTION 

(AN X GRID LINE EVERY 4th POINT 
(HAKES GRID UNIFORM 



IRFADS NO-PARTICtE DATA 
• NOT AN OVERLAPPING FRANSFFR 
i Xyz SETS DISK FILE POINTFR TO 
(EITHER MOTOR OR EXHAUST DATA 
(READS PARTICLE DATA 



(JUST GOOD PRACTICE TO CIOSE 
(I/O PATHS 



5040 MOVE D(I,2) ,0 

5050 LABEL USING *.DDDD*;D(I,2) 

5060 NEXT I 
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! NUMBER Y AXIS 




50*50 FOR 1=10 TO Max STEP 10 

5100 MOVE D(4,2)-G,I 

5110 LAPEL USING , *,K’jI 

5120 NEXT I 

5130 CSIZE 6,. 6 

5140 Titlp»= , TUO-AUGLE METHOD- 

5150 Title1t='For Various Angle Ratios* ! STRINGS FOR PI OTS 

51 fcO Si7et=‘SlZE (nicrons)' 

5170 Thtat="THETA (rad) 1 

5130 Siib$=T 

5170 LOIR 90 

5200 LORG 5 

5210 B=D(4,2)-G-(D(E,2)-D(1,2))/10 

5220 MOVE R, Max/2 

5230 LABEL Size* 

5241 LOIR 0 

5250 A=(D(E,2)+D(1 ,2))/2 

5260 HOVE A, -Max/4 

5270 LABEL Thtai 

5280 LORG 3 

5290 HOVE A, -Max/4 

5200 LABEL Sibt 

5310 LORG 5 

5320 MOVE A,Maxtl.l5 

5330 LAPEL Tit let 

5340 CSIZE 4.5 

5350 MOVE A,Max«l .05 

5360 LABEL Titlelt 

5370 PEN1JP 

5380 SUP-END 

5390 SUB BistribmionlDU) ,Tratio,E) 



5400 


LINE TYPE 1 




5410 


CLIP ON 




5420 


Tratiol=VALS(Tratio) 


! ANGLE RATIO 


5430 


MOVE 0(1 ,2), 0(1,1) 




5440 


FOR 1=1 TO E 


! THIS SUBROUTINE PLOTS THE 


5450 


DRAW D(I,2) ,D(I,1 ) 


'PARTICLE SIZES DERIVED USING 


5460 


NEXT I 


! VARIOUS ANGLE RATIOS APPLIED 


5470 


LORG 4 


! TO THE DATA OVER A RANGE OF 


5430 


LINE TYPE 1 


(ANGLES 


54?0 


CSIZE 4 ,.3 




5500 


HOVE D(I-1 ,2) ,D(I-1 ,1 )+2 




5510 


LABEL Tratiot 




5520 


PENUP 





553) SUBEND 
5540 SUB Tuoangle 
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5550 

5560 

5570 

5590 

5590 

5600 

5610 

5620 

5630 

5640 

5650 

5660 

5670 

5680 

5690 

5700 

5710 

5720 

5730 

5740 

5750 

5760 

5770 

5780 

5770 

5900 

5B10 

5«2H 

5830 

5840 

5850 

5860 

5870 

5880 

5890 

5900 

5910 

5920 

5930 

5940 

5950 

5960 

5970 

5990 

5990 

6000 

6010 

6020 



ix***xxx*****x****x***#******xx****xxxxxxx*xx*xx*x*xx*xx*x<**i*xxx 
ixxxxxxxK'mmx SlIPPPOGRAM ' TWOANGLE' xxxixxxxxxxxxxxxxxx 



(xxxxxxxxxxxxxxx 
ixxxxxxxxxxxxxxx 
ixxxxxxxxxxxxxxx 
! xxxxxxxxxxxxxxx 



XXXXXXXXXXXXXXXXXXX 

XXXXXXXXIXXXXXXXXXX 

XXXXXXXXXXXXXXXXXXX 

XXXXXXXXXXXXXXXXXXX 



PARTICLE SIZING BY 
MEASURING THE RATIO 
OF INTENSITY AT 
TWO ANGLES 

(xxxxxx xxxxxxxxxxxxxxx* ****** xxxxxxxxxxxxxxx***** xxxxxxxxxxxxxxx** 

OPTION BASE 1 

COM /Tup / AvlUI)?4),H,h1 ,F 

CON /Gauss/ T1(1024),L ! THE GAUSS IS NOT USED HERE 

• PUT THE CON BLOCK HAS THETA 

DIM D(200,2) ! AND WAVELENGTH 

PRINT CHR* (12) 

PRINT USING •///////* 

PRINT ’ THIS SUBPROGRAM USES. THE TWO-ANGLE METHOD DESCRIBED BY BUCHFIE* 
PRINT ■ ‘ 

PRINT * TO CALCULATE PATICLE SIZE FOR VARIOUS ANGLE RATIOS AND ANGLES,* 
PRINT * * 

FRINT * IT IS HOPED THAT THE CURVES WHICH RESULT WILL SHOW * 

PRINT * 8 

PRINT * WHICH SIZE IS THE MOST PROBABLE.* 

PRINT * " 

PRINT ‘AFTER NOTING FROM THE RAW DATA WHICH ANGLES CONTAIN THE CFNTER LOBF- 
PRINT *“ 

PRINT * ENTFR»************THE SMALLEST USEABLE ANGLE,* 

PRINT * THE SMALLEST ANGLE RATIO, AND" 

PRINT " THE STEP SIZE BETWEEN ANGLE RAIOS* 

PRINT * YOU WISH TO EXPLORE* 

PRINT * EXAMPLE************* ,012, 1.2,. 4* 

INPUT * ENTER thetal, theta-ratio, theta-step*, Q, A. B 

OFF KEY 

X=0 (GRAPH WILL APPEAR ON SCREEN 

Begin; '.RUNNING CONTINUES HERE WHEN A HARD COPY IS DESIRED 

! Van Dp Hnlst and Gvnprech 6 Sleppevich explain that the chanoe in 
(wavelength of the bean is accounted for by dividing by the index 
! of refraction of the nediun. x***xx***L=L/Hl xxxx****** 

C=f L/M1 / . 5?/P I ) *2 ! see page 15 of nasa tech oaoer 2156 

FOR N=1 TO 1024 !to see this is a convenient constant 



IF QUKN) THEN 5950 
NEXT N 

FOR Tratio=A TO 3 STEP B 

FOR J=N TO 1010/Tratio STEP 10 

ThtaWKJ) 

Thta2-Tl (J*Tratio) 
Deltatheta=Thta2‘2-ThtaP2 
11*0 
12=0 



! FINDS POSITION OF HINIHUN ANLGE 

(VARIOUS ANGLE RATIOS 
(SETS THE RANGE OF POSSIBLE 
! ANGLES (THETA 1) 

! THETA 2 
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6050 

6060 

6070 

6080 

6090 



T2=I2+Avl(I) 
NEXT I 

I2-T2/11 

E=(J-N)/10+1 



FOR I=INT(J*Tratio-5) TO INT(J*Tratio+5> 



NEXT I 

11 = 11/11 



i DONE HERE FOR THETA2 FOR 
! THE GIVEN ANGLE RATIO 



! OF A STEADIER CALCULATION 



6100 



6110 

6120 

6130 

6140 

6150 

6160 

6170 

6180 



IF 1 1 < 12 OR 11=0 OR 121=0 THEN D(E,1)=0 ! ALLOWANCE FOR IF THE DATA 
IF II <12 OR 11=0 OR I2<=0 THEN GOTO Xconp IIS NOT WELL BEHAVED 



D(E, 1)=S9R( -C/Del tatheta *LOG{I2/I1 )) 



■THIS IS A COUNTER FOR THF ARRAY CONTAINING 



'PARTICLE SIZE ANO THETA1 FOR THE GIVEN 
'ANGLE RATIO 



! VALUE OF DIAMETER based on INTENSITY RATIO 
IF or a given ANGLE RATIO and ANGLE THETA 1 



I TWO ANGIE METHOD 



6190 Xconp: 



B(E,2'=Tbt?1 

Spar=PI*D(E,l)*Hl/L (THIS IS 
Tbar=SDar*lbta2 
IF Tbar >4 THFN J=1010/Tratio 
IF E=1 AND Tbar >3 THEN 6330 



6200 

6210 

6220 

6230 



PI*D/LAMBDA - THE SIZE PARAMETER 
I THETA BAR F OR THF LARGER ANGLE 
I THIS ENDS THE DO LOOP FOR THIS 
(ANGLE RATIO SINCE THE GAUSSIAN 



6240 (the above line ends all calculation if IIS NOT VALID WHEN Tbar > 3 
6250 'the first element (E=l)failed the test 
6260 NEXT J 



6290 PRINT CHPFI12) 

6300 IF Tratio=A TH r N CALL Plot2(E,D(»),X> 

6310 CALL Distributional*) ,Trat io,E) 

6320 NEXT Tratio 

6330 ON KEY 2 LAPFL "HARD COPY" GOTO Hard 

6340 ON KEY 3 LABEL "HENUE ’ GOTO Subexit 

6350 PRINT USING •/////■ 

6360 PRINT * YOU CAN GET A HARD COPY BY PRFSSING KEY t 2" 

6370 PRINT ' OR EXIT THIS ROUTINE BY PRESSING KEY » 3" 

6380 Standbv : GOTO Standby 

6390 Hard : X=1 !A VARIABLE TO CONTROL PLOTTER 

6400 GOTO Begin 

6410 Subexit: GINIT 

6420 GRAPHICS OFF 

6430 SUBEND 

6440 SUB 9iift 

6450 COM /Readata/ B,P,H,Q3$[20I,Q4$[201,Zz$t20],Y1<») BUFFER, Y2<*> BUFFER 
6460 DIM E< 1 : 8 1 92) 

6470 PRINT CHR$( 12) 

6480 PRINT USING "///////' 

6490 PRINT " RAW DATA IS BEING SHIFTED TO CORRECT FOR SHALL GAPS BETWEFN SCANS" 
6500 PRINT ’* 

6510 PRINT " BE PATIENT, IT'S A LONG SET OF DO-LOOPS" 



6270 

6280 



FIRST TIMF THROUGH- D HAS THE MOST ELEMENTS IT WILL HAVE 
SO SET UP THE GRAPH USING D's PRESENT PARAMETERS 
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6520 

6530 

6540 

6550 

6560 

6571 

6580 

6590 

6600 

6610 

6620 

6630 

6640 

6650 

6660 

6670 

6680 

6690 

6700 

6710 

6720 

6730 

6740 

6750 

6760 

6770 

6780 

6790 

6800 

6810 

6920 

6830 

6840 

6850 

6860 

6870 

6880 

6890 

6900 

6910 

6920 

6930 

6940 

6950 

6960 

6970 

6980 

6990 

7000 



THFRE ARE SOME SHALL GAPS BETWEEN SCANS AND THIS SUBROUTINE 
SHIFTS THE DATA SO THAT THE FIRST DIODE DATA POIMt IS HOvED 

TO THE VERY BEGINNING OF ITS 1024 BLOCK IN THE OVERALL ARRAY 
. THE FIRST SET IS RIGHT ON, THE NEXT IS ONE OFF, THE THIRD IS 
TWO OFF ,,SO FORTH. THIS HAY NOT HATTER WITH OUR RESOLUTION 
AND IS PROBABLY DIE TO THE HEHORY CARD CYCLING AT THE END OF 
SOME SCANS. SEE THE CIRCUIT TIMING DIAGRAM IN THE THESIS. 



SELECT B 
CASE 4 
H=3 

CASE 8 

Y\-7 

END SELECT 
FOR KM TO 2 



! EXHAUST DATA HAS 4 SCANS 



! MOTOR DATA HAS 8 SCANS 



IF KM THEN HAT E= Y1 ! OPPFRATES ON NO-PARTICLE AND PARTICLE SETS 

IF K=2 THEN HAT E= Y2 
FOR J=0 TO H 



IF H=7 AND J<=3 THEN 6800 'One 4096 Block Doesn't Need Shifting 
FOR I=< J ) »1 024+1 TO (J+t 1*1024 ! BLOCKS OF 1024 

IF H=3 THEN L=I+J lArray B Has the worst 

IProblen With Shifting Data 

IF H=7 THEN LM+1 ! Array D is Allways off bo one 
IF LMJ+DM024 THEN LMJM)*1024 
•JUST TO AVOID PROGRAHHING ERROR AT THE END OF THE APR AY 
E(I)=E(L) ITHIS SHIFTS THE DATA 

• DEPENDING OH 'l ARRIVED AT 

• BY OBSERVING RAW DATA 

NEXT 1 
NEXT J 

IF KM THEN HAT Yl= E 
IF K=2 THEN HAT Y2= E 
NEXT K 
SUBEND 

SUB Store ! THE LARGE ARRAY CONTAINING MULTIPLE SCANS HAS BEEN REDUCED 
i TO A HEAN SCATTERING PROFILE BY AVERAGING AND FILTERING. IF 
! YOU FEFL THAT STORING THIS DATA IS NECESSARY THIS ROUTINE 
i DOES IT. II SAVES DISK SPACE TO STORE THE RESUI IS THEN 
! El IH1NATE THE RAW DATA IF YOU FEEL CONFIDFNT THAT THE 
i REDUCTION IS THE BEST THAT CAN PE. IN OTHER WORDS, 

! D" NOT PURGE A RAW DATA FILE UNLESS YOU ARE ABSOLUTELY SURE 
! YOU WON'T WANT TO REDUCE IT AGAIN. 

COM /Hrdqans/ Ao2(») !Av2(») IS THE REDUCED DATA 

DIH DataU :1024) BUFFER ! BUFFER USED TO TRANSFER TO DISK 

HAT Data= Av2 
PRINT CHR t ( 1 2 ) 

PRINT USING *//////* 

PRINT ’ A SUGGESTED HETHOD OF NAHING REDUCFD DATA FILES IS AS FOLLOWS" 

PRINT * H HOT OR BE AH* 

PRINT * X EXHAUST BEAH" 
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7010 

7030 

7030 

7040 

7050 

7060 

7070 

7080 

7090 

7100 

7110 

7120 

7130 

7140 

7150 

7160 

7170 

7180 

7190 

7200 

7210 

7220 

7230 

7240 

7250 

7260 

7270 

7280 

7290 

7300 

7310 

7320 

7330 

7340 

7350 

7360 

7370 

7380 

7390 

7400 

7410 

7420 

7430 

7440 

7450 

7460 

7470 

7483 

7490 



PRINT * C CALIBRATION , IF NO 'C' THEN AN ACTUAL FIRING IS ASSUMED* 

PRINT " G GLASS BEADS 1 

PRINT * IF Nn 'C' THEN 'G' STANDS FOR 'GAP' PROPELLANT* 

PRINT * A ALUMINUM OXIDE* 

PRINT * PPP CHAMBER PRESSURE FOR RUN OR OTHER. ... (NOZZLE TYPE)" 

PRINT * PARTICLE SIZE FOR CALIBRATION* 

PRINT * MHDD MONTH, DAY OF RUN OR CALIBRATION* 

PRINT ** 

PRINT * EXAMPLE: HCA125JN12* 

PRINT ** 

PRINT * MOTOR BEAM CALIBRATION USING ALUMINUM OXIDE FROM 1 TO 25 MICRONS* 
PRINT * ON JUNE 12" 

PRINT * EXAMPLE: XG550AU10* 

PRINT *" 

PRINT * EXHAUST BEAM DATA FOR GAP PROPELLANT AT 550 psi ON AUGUST 10" 
PRINT " PLACE A DISK IN THE RIGHT HAND DRIVE" 

PRINT " ENTER THE FILENAME YOU NISH TO USE FOR THIS REDUCED DATA* 

INPUT At 



CREATE BDAT At, 512,16 

ASSIGN PDisk TO At 

ASSIGN Muff TO BUFFER DataU) 

CONTROL Muff, 3;1 ,8192,1 '.BUFFER IS FULL 

TRANSFER ?Suff TO gDisk jCOUNT 8192 
WAIT FOR FOT BDisk 
ASSIGN Muff TO * 

ASSIGN RDisk TO » 

SUPEND 

SUB Review(AvlU)) 

DIM DataU: 1024) PUFFER 

PRINT * EACH REDUCED FILE CONTAINS ONE SCAN ONLY . YOU MUST KNOW IF IT* 
PRINT * IS EXHAUST , MOTOR CAVITY , OR CALIBRATION DATA.* 

PRINT ** 

PRINT * THE DISK WITH THE REDUCED FILE SHOULD BE IN THE RIGHT-HAND DRIVE" 
PRINT *" 

PRINT * ENTER THE FILENAME OF THF REDUCED DATA" 

PPINT " TO BE REVIEWED* 

PRINT *" 

1NFUT At 

ASSIGN PDisk TO At 

ASSIGN Muff 10 BUFFER DataU) 

CONTROL SDisk,5;l 

CONTROL Muff, 3;1, 0,1 ITHIS IS AN EMPTY BUFFER 

TRANSFER 8 D 1 U TO MuffjCOUNT 8192 ! 1024*8=NUM9ER OF BYTES 

WAIT FOR EOT POisk 
ASSIGN PDisk TO * 

ASSIGN Muff TO » 

HAT Avl= Data 
SUBEND 

SUB Fi 1 ter ( A < *) ,Fi 1 ) 
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wum. 




7520 


PRINT USING •////////' 




7530 


PRINT " 


FILTERING THE SCATTERING PROFILE 1 


7540 


FOR M TO Fil 




7550 


FOR 1=1 TO 1014 




7560 


B=I+1 


! THIS IS A Synetric Moving Average 


7570 


C=I+2 


! TYPE OF DIGITAL FILTER. EACH 


7580 


8=1+3 


'.DATA POINT IS EQUALLY WEIGHTED 


7590 


F=I+4 


•IN THIS CASE PUT THIS CAN PE 


7600 


G=I+5 


! CHANGED IF ONE DETERMINES THAT 


7610 


H=I +6 


(FEWER POINTS WITH UNEQUAL WEIGHTS 


7620 


K=I+7 


IUOULD BE FASTER OP GIVF FETTER 


7630 


L=I+8 


■RESULTS. THIS TYPE OF FILTER WAS 


7640 


N=I+9 


(USED SINCE IT INTRODUCES NO PHASE 


7650 


P=IM0 


'LAG (Angular Resolution) . 


7660 


E(Q)=<A(I)+A(B)+A(C)+A(D)+ACF)+A<G)+A(H)+A(K)+A(L)+A(N)+A(P>)/11 


7670 


NEXT I 




7680 


HAT A= F 




7690 


NEXT J 




7710 


SUPEND 
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